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Automatic contrast medium control in images 



5 

The present invention relates to the field of digital imaging. In particular, 
the present invention relates to a method of processing a series of x-ray images of an 
object of interest, to an image processing device and to a computer program for 
processing a series of x-ray images of an object of interest. 

10 In diagnostics and therapy of, for example, a coronary heart disease, 

projection x-ray images acquired in so-called cine imaging mode play a crucial role. To 
make the coronary arteries visible in such coronary angiographies, a radio-opaque 
contrast medium is applied by means of a catheter placed, for example in the entrance to 
a main coronary artery. Pre-interventional coronary angiographic sequences showing 

1 5 the coronary vessel tree during several heart cycles serve as diagnostic images to detect 
stenoses, and as roadmaps for the intervention itself. 

During the intervention, a catheter or a guide wire is advanced under x- 
ray monitoring through the vessel to the lesion. Only occasional bursts of contrast agent 
for verification purposes can be given while this procedure is performed. To help 

20 navigation, a single frame showing the entire vessel tree filled with contrast agent is 
therefore selected manually from the pre-interventional angiograms to serve as a 
roadmap and is displayed on a screen next to the interventional live images. This 
roadmap, however, is naturally static and is hence not consistent with the instantaneous 
heart and respiration status in the live images. Due to this, a navigation of the catheter 

25 or the guide wire is difficult since the roadmap does not match the interventional image. 

It is an object of the present invention to provide for an improved 
processing of a series of x-ray images of an object of interest, where the object of 
interest is visible due to a contrast medium. 

According to an exemplary embodiment of the present invention as set 

30 forth in claim 1, the above object may be solved with a method of processing a series of 
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x-ray images of an object of interest, wherein the object of interest is visible due to a 
contrast medium, wherein, according to the method, an image of the series of x-ray 
images where the object of interest is not sufficiently filled with the contrast agent is 
automatically determined. 
5 Due to this automatic determination, images, for example, of the 

coronary angiograms can be determined where the vessel tree is sufficiently filled with 
a contrast agent. Such images may then be preferably displayed as live images, for 
example, for catheter navigation. However, a couple of images where the vessel tree is 
sufficiently filled with contrast agent may be overlaid to produce an image where the 

10 complete vessel tree is visible. Also, images may be automatically determined where the 
complete vessel tree is filled with contrast agent. 

According to another exemplary embodiment of the present invention as 
set forth in claim 2, a pre-processing of the images is performed such that a background 
of the object of interest is at least partially suppressed. Furthermore, enhancement of 

1 5 part of the object of interest may be performed. 

According to another exemplary embodiment of the present invention as 
set forth in claim 3, a morphological filtering is performed and an accentuation of part 
of the object of interest visible in the respective images of the series of x-ray images is 
performed on the basis of a motion of the object of interest. Preferably, this may allow 

20 for an improved image quality. 

According to another exemplary embodiment of the present invention as 
set forth in claim 4, first and second order derivatives of images are used to enhance 
image information relating to the object of interest. 

Furthermore, this may allow to further improve an image quality such 

25 that, for example, in coronary angiograms the vessel tree displayed has an improved 
image quality. 

According to another exemplary embodiment of the present invention as 
set forth in claim 5, a determination with respect to whether the object of interest is 
sufficiently filled with contrast agent is performed on the basis of a number of picture 
30 elements of the image having a pixel value exceeding a pre-set threshold value. 
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Advantageously, this allows for a fast and robust determination of 
images, where the object of interest is sufficiently filled with contrast medium. 

According to another exemplary embodiment of the present invention as 
set forth in claim 6, a determination of the image of the series of x-ray images, where 
5 the object of interest is not sufficiently filled with the contrast agent is performed on the 
basis of a histogram analysis, a feature curve analysis and a feature curve segmentation. 
The feature curve segmentation may be based on a maximum likelihood segmentation. 

This may allow for a fast, efficient and accurate determination. 

According to another exemplary embodiment of the present invention, a 
10 95-percentile of a histogram is used to determine a feature function. Then, a second 

histogram is formed on the basis of the feature function. Then, a threshold is determined 
in the second histogram, which separates a first state of an image where the object of 
interest is insufficiently filled with contrast medium, from a second state of an image 
where the object of interest is sufficiently filled with the contrast medium. 
1 5 According to another exemplary embodiment of the present invention, 

the method is adapted for determining images of coronary angiograms where the vessel 
tree of the heart is sufficiently filled with contrast agent. 

Claim 1 1 relates to an image-processing device according to another 
exemplary embodiment of the present invention, which may allow for an improved 
20 operation, for example, in conjunction with coronary angiograms and/or catheter 
navigation. 

The present invention relates also to a computer program, for example, to 
an image-processing device for processing a series of x-ray images of an object of 
interest. The computer program according to the present invention is defined in claim 

25 12. The computer program according to the present invention is preferably loaded into a 
working memory of a data processor. The data processor is thus equipped to carry out 
the method of the invention. The computer program may be stored on a computer 
program medium such as a CD-ROM. The computer program may also be presented 
over a network such as the Worldwide Web, and can be downloaded into the working 

30 memory of a data processor from such a network. The computer program may be 
written in any suitable programming language, such as C++. 
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In the following, what might be seen as the gist of an exemplary 
embodiment of the present invention is described with respect to the exemplary 
identification of complete vessel tree in coronary angiograms. First, a feature map 
where the vessel regions have been enhanced is created. Then, a vessel enhancement is 
5 based on the above observations and vessels are the locally darkest oriented structures 
with significant motion. A background is then equalized by a morphological top hat 
transformation with a structuring element comparable to the diameter of the last vessel. 
Then, a vessel motion is exploited to increase the contrast between background and 
vessel. Vessel borders are amplified by calculating a gradient norm, while vessel centers 
10 - which do not respond to first derivatives - are enhanced by a second derivative 
operator. The resulting enhanced images can be regarded as containing two classes, 
namely bright vessels and dark background. It may therefore also serve as a basis for 
vessel segmentation. 

According to an exemplary embodiment of the present invention, a 
1 5 histogram of the enhanced images is determined and the 95-percentile of these 

histograms is used as a feature. The more area in an image is covered by contrast agent, 
the higher the 95-percentile. When plotted over time (i.e. the frame number of the 
images), the feature curve clearly shows the following phases: an in-flow of the contrast 
agent; a filled state, where the vessel tree is sufficiently filled with contrast agent, and 
20 an out-flow, where the contrast agent is washed out of the vessel tree. This feature curve 
is then segmented. 

According to an exemplary embodiment of the present invention, the 
observations in each of these three states is modeled by a polynomial and then, a 
segmentation is performed, which allows the best fit of three polynomials as measured 
25 by a maximum likelihood criterion. 

According to another exemplary embodiment of the present invention, 
the state sequence is modeled by a Hidden Markov model and is estimated by a 
maximum a-posteriori (MAP) criterion. 

These and other aspects of the present invention will become apparent 
30 from and elucidated with reference to the embodiments described hereinafter. 
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Exemplary embodiments of the present invention will be described in the 
following, with reference to the following drawings: 



Fig. 1 shows a schematic representation of an image-processing device 
5 according to an exemplary embodiment of the present invention, 

adapted to execute a method according to an exemplary 
embodiment of the present invention. 
Fig. 2 shows a block diagram of a spatio-temporal filtering chain 

according to an exemplary embodiment of the present invention. 
10 Fig. 3 shows a processing chain with first and second order derivatives 

as implemented according to an exemplary embodiment of the 
present invention. 

Fig. 4 shows vessel map logarithmic-histograms for an in-flow frame (a) 
and a filled state (b) according to exemplary embodiments of the 
1 5 present invention . 

Fig. 5 shows a percentile feature of a frame index (a) filtered by a FIR 
low pass filter of length 7 (b) according to an exemplary 
embodiment of the present invention. 
Fig. 6 shows a histogram of the 95-percentile feature curve depicted in 
20 Fig. 5b according to an exemplary embodiment of the present 

invention. 



In the following, the present invention will be described with respect to a 
statistical model based identification of complete vessel tree frames in coronary 
25 angiograms. However, it should be noted that the present invention is not limited to this 
application, but may be applied to any imaging applications where the object of interest 
is made visible by application of a contrast medium such as x-ray imaging, ultra sonic 
imaging or MR imaging. 

Fig. 1 depicts an exemplary embodiment of an image processing device 
30 according to the present invention, for example, executing an exemplary embodiment of 
the method in accordance with the present invention. The image-processing device 
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depicted in Fig. 1 comprises a central processing unit (CPU) or image processor 1 
connected to a memory 2 for storing time series x-ray images. The image processor 1 
may be connected to a plurality of input/output network or diagnosis devices, such as an 
x-ray scanner. The image processor is furthermore connected to a display device 4 (for 
5 example, a computer monitor) for displaying information or images computed or 

adapted in the image processor 1 . An operator may interact with the image processor 1 
via a keyboard 5 and/or other input/output devices, which are not depicted in Fig. 1. 

According to an exemplary embodiment of a method of operating the 
image processing device depicted in Fig. 1, the image processing device firstly performs 
10 a vessel enhancement. The vessel enhancement firstly comprises a morphological 
filtering and then an exploitation of the motion of the vessel tree. Then, vessel 
information is enhanced in the angiograms by using first and second order derivatives of 
the angiograms. 

Then, after the vessel enhancement, a histogram analysis and feature 
15 curve analysis is performed. 

In a next step, feature curve segmentation is performed. The feature 
curve segmentation may perform maximum likelihood segmentation and/or 
segmentation on the basis of a Hidden Markov model. 

This will be described in further detail in the following: 
20 Vessel enhancement 

Morphological filtering 

A first step towards enhancing vessels with respect to background 
exploits the property that, being filled by a radio-opaque contrast agent, vessels are 
locally darker than their immediate surroundings. In order to flatten the varying 

25 background intensity, a background image is calculated by removing foreground 

structures. To this end, a local sliding maximum filter is firstly applied, the dimensions 
of which are slightly larger than the largest expected vessel diameter. This operation 
removes locally dark structures, which are smaller than the size of the filter. However, 
edges of larger dark regions bordering on brighter regions are also eroded, i.e. moved 

30 towards the inside of the dark region by approximately the size of the filter. An example 
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of such a region in coronary angiograms is the diaphragm, which usually appears as a 
dark half disc in the lower portion of the angiograms. 

To reconstruct such eroded edges, a local sliding minimum filter of the 
same size may be applied, which propagates the boundary back to its approximate 
5 original location. Maximum and minimum filters may be applied in windows of a pre- 
determined size and implemented as successive horizontal and vertical ID filters in 
order to save complications. 

The filtered images then contain almost no foreground information. 
Subtraction of the original frame from its maximum - and minimum - filtered 
10 background version thus leaves the foreground vessel information. The succession of 
maximum and minimum filtering may be referred to as morphological closing and 
taking the difference between the original and its closing may be called top hat filter. 
Due to the maximum filter applied in the first filtering step, the gray level at each pixel 
in the closing-filtered image is never smaller than the gray level of the same pixel in the 
15 original. The gray levels in the top hat filtered image are therefore non-positive. 

Exploiting motion 

As has been found out according to the present invention, an additional 
clue to vessels is their mostly strong and jerky motion, which is caused by the beating 
heart. When the vessel moves to a new position, absorption at this position will increase 

20 due to the contrast agent within the vessel. Therefore, when calculating the pixel-wise 
difference image between any given frame and its predecessor (or between the current 
frame and a frame some fixed distance back in time), pixels with the new vessel 
positions will tend to exhibit negative differences, while the vessel positions in the 
previous frame will tend to show positive differences. This effect is even more 

25 pronounced in top hat filtered angiograms, since the top hat filter reduces background 
structures, which would cause clutter in the difference images when moving themselves. 
Therefore, positive values in the difference image are clipped to zero and then, the 
clipped difference image is added to the current top hat-filtered frame. The overall result 
is a tendency to make the vessels even darker. Local minima of small extent, which is 

30 unlikely to be caused by moving vessels, can optionally be removed by a morphological 
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closing, for example, with a 3 x 3 structuring element. A block diagram of the 
processing chain described above is shown in Fig. 2. 

The block diagram of the spatio-temporal filtering chain described above 
is depicted in Fig. 2. X n denotes an n-th incoming angiography frame and Y n its top hat 
5 filtered version wherein the top hat filtering is performed in top hat filter 24. After 

subtracting its predecessor Y n -i (Y n _i is hereby transmitted from z'-element 25) from Y n 
by means of subtracting element 21, positive differences are clipped towards zero in 
clipping element 26. After optional weighting by means of weighting element 22 and 
closing by means of a 3x3 closing element 27, the result is added to Y n by means of 
1 0 adding element 23 . 

First and second-order derivatives 

To enhance vessel information in the angiograms after spatio-temporal 
processing by top hat filtering and motion analysis, a combination of first and second 
order derivatives may be used to enhance vessel information. Taking a magnitude of the 

15 first derivative of the angiograms images evidently captures the steep slopes at the 

vessel borders, but naturally fails to respond to the flatter portions of the profile within 
the vessel. Therefore, the first order derivative is complemented by a second-order 
derivative. Since the second-order derivative captures the middle of the vessel, both 
derivatives are added to obtain a response over the entire vessel profile. 

20 In a 2D-image f(x 9 y), the absolute first derivative is replaced by the 

norm |A f(x, y)| of the gradient. Gradient information is computed from the eigenvalues 
XI and X2 of the 2 by 2-tensor T introduced for orientation analysis. Apart from the 
gradient, the eigensystem of this tensor provides also information on the local structure 
of the image signal, like orientation, homogeneity and corners. In continuous 

25 coordinates, the tensor T is given by 



T= f (V/)(V/) r rffi (1) 
Ja 

where £1 is a small local window. In discrete coordinates, the horizontal 
and vertical derivatives are first calculated using symmetric finite difference kernels. 
30 For the entries along the main diagonal of T, these derivatives are squared, while the 
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entries on the counter-diagonal are given by the pixel-wise product of horizontal and 
vertical derivatives. Finally, the integration over Q is realized by a lowpass filter of 
appropriate size. The gradient norm is then given by the square root of the sum of the 
eigenvalues, i.e. 

5 

I V/(* f y)\ = y/*i{**v) + **fav) (2) 

In addition, the eigenvalues can be used to distinguish three different 
types of structure: if both eigenvalues are (close to) zero, the local structure in the 
support window Q. is flat. If the eigenvalue XI is large, while the smaller eigenvalue A2 

10 is (close to) zero, i.e. rank(T) = 1, there is a structure with a unique orientation in Q, 
caused, e.g., by a locally straight vessel. In this case, the eigenvectors run parallel 
respectively perpendicular to the orientation. Finally, if both eigenvalues are relatively 
large, T is not rank-deficient. Then, there is only one local orientation, caused by a 
sharply bending vessel projection. The gradient norm fails to respond to the vessel 

15 middle. Therefore, the gradient norm is complemented by the Laplacian 



A/Ovy) = fasffc v) + 5jff y) = M x * y) + fni x > v) 

which is linear, shift-invariant and, in its continuous version, isotropic. 
This operator hence responds to vessel ridges independent of their orientations. As a 

20 discrete approximation to the Laplacian operator, a difference of Gaussians (DoG) filter 
is used, which is an approximation to a combination of the Laplacian with a Gaussian 
lowpass filter, or so-called Marr-Hildreth operator. The result of filtering the top hat 
prefiltered image by this operator is added pixel-wise to the gradient-norm image. 

Both computation of the gradient norm with a calculation unit 33 via the 

25 tensor T and the second derivative with another calculation unit 34 involve lowpass 

filtering, which is e.g. performed by lowpass filter 35. On the one hand, this makes the 
processing results more robust with respect to noise, while on the other hand, it blurs 
vessel edges. The top hat filter, however, preserves edges. In a final step, the sum 
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image, which is provided by adding element 31, is therefore multiplied by the top hat- 
filtered image by means of multiplier 32 to reduce the blur of the vessel edges. A block 
diagram of this processing chain is depicted in Fig. 3. 

Histogram analysis and feature curve 
5 As found according to the present invention, a presence of contrast agent 

increases a frequency of occurrence of bright intensity in the vessel map. Therefore, 
vessel map histograms are analyzed over the time - i.e. frame index - to obtain a time 
dependent feature curve related to the presence of contrast agent. Since vessels filled by 
contrast agent cover about 5% of the total image area, as has been found according to 

10 the present invention, a 95-percentile of the histograms was chosen as scalar feature. 

The presence of contrast agent increases tails of the vessel map histograms and thus also 
the 95-percentile. The three states in-flow, filled state and out-flow are thus visually 
readily identifiable in the feature function formed over the frame index. 

This may be taken from Figs. 4a and 4b and from Figs. 5a and 5b. Figs. 

15 4a and 4b show logarithmic vessel map histograms for an in-flow frame (a) and a filled- 
state frame (b). Thus, two different states are depicted. It should be noted that the 
present invention is also applicable for a plurality of states, for example, in a sequence. 
E.g. an inflow state, a filled state and an outflow state may be determined. The 95- 
percentile in Figs. 4a, 4b, 5a and 5b is given by the arrows. 

20 As may be taken from Figs. 4a and 4b, the 95-percentile in Fig. 4a is 

approximately 45 and the 95-percentile in Fig. 4b is approximately 63. The abscissae in 
Figs. 4a and 4b map the gray values in the image of interest. The ordinates of Figs. 4a 
and 4b map the frequency of these gray values. The abscissae in Figs. 5a and 5b relate 
to the frame numbers of the images and the ordinates relate to the 95% percentile 

25 determined on the basis of a histogram of the preprocessed image as in Figs. 4a and 4b. 

The evolvement of the percentile feature over time is depicted in Fig. 5, 
where Fig. 5a depicts the percentile feature over the frame index and Fig. 5b shows the 
same percentile feature over the frame index but filtered by a FIR low pass filter of 
length 7. 

30 As may be taken from Figs. 4a and 4b and Figs. 5a and 5b, two states, in- 

flow and filled state, are visually readily identifiable. Also, it may be taken from the 
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Figs. 4a, 4b, 5a and 5b that the area covered by the vessels - and thus the percentile - 
depends on the heart cycle, since the filled state exhibits a periodic oscillation. 
Feature curve segmentations 

Since each of the states in-flow, filled state and out-flow must form a 
5 coherent region, the feature curve as depicted in Figs. 5a and 5b, but preferably the 
filtered feature curve depicted in Fig. 5b, is segmented into three regions. For a feature 
curve of length N, there exists 



M = (2j = (N-2)!2? "-r & 

10 different segmentations. For a typical length of N = 70, one obtains N = 

2415 possible segmentations. 

Maximum likelihood segmentation 

One possibility to segment the feature curve is by least-squares fitting 
three polynomials of a given degree to the feature curve. Under Gaussian assumptions, 
15 this is equivalent to a Maximum-Likelihood (ML) segmentation. The use of ML- 
estimation rather than seeking the Maximum a Posteriori (MAP) estimate may be 
justified by the above strong coherency constraint on the solution. Denoting the 
succession of states by Q = {q } , 9 q N }, where q* is the state label for the fth frame of 

the angiogram, and the sequence of percentiles by S ={y,,...,^}, with Si being the 
20 percentile observed for the ith angiogram frame, the solution Q which obeys 



3 

Q = arg maxp(5|Q) = org max |"| p(s } \mj. aj) (5) 
is sought. 

Here, j - 1, 2, 3 is the index of region Rj , and Sj is a vector with ordered 
25 observations in each region. Modeling the observations as Gaussian distributed, mj and 
a j are the region-specific parameters of the distribution. Also, the realistic assumption 
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has been made that the observations are region-wise independent. Furthermore, it is 
assumed that the observations within each region are conditionally independent, i.e. that 
their dependencies are fully captured by the coherency of the region. Then, for the 
region likelihood, the following may be obtained: 

Estimating the unknown region parameters mj and a] by 



10 where Nj is the number of frames in region Rj , and replacing the true 

parameters by their estimates, Eq. (6) yields 



pCsylt&j.aJ) = y/^d* ' -e"* (8) 
Thus, the following equation is used 

15 

= a rg mm ^ ■ log (&f) (9) 

where all constants and factors not influencing the maximization have 
been dropped, and where it has been have taken into account that Nl +N2 +N3 = N. 
Evidently, this formulation seeks the segmentation which permits fitting a polynomial 
20 of degree zero optimally - i.e. with minimum mean square error as defined by Eq. (7) - 
to each region. 

An alternative is to segment such that a polynomial of degree one can be 
fit best to each region. Then, each region mean mj is replaced by the slope aj and the 
mean bj , and the likelihood for each region becomes 



:aigoiux J]l<Sj|*fcj.^ 



arc* max 



n 
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pfo|«,,fr, )C r?) = ( . 1 J exp{--i T (*•* - «j77. - tj) 2 } 



(10) 



With the estimates 



* 2 



Tit 

one again gets 



pfola^.oj) ~ y^Jf (12) 
The parameters <5 y and are estimated by minimizing <3"j : 



0 =s> *» = a * n + N ih ( 13 ) 



= 0 a* n ^ =i J X] « 3 + &j « (14) 

n^Itj n£Jlj tiGJIj 



These equations, which must be solved for a y and b } , can be decoupled 
when centering the region-internal coordinate n such that 

X> = 0 (15) 
15 This yields 



^-l^E*" «* = £ ' » 2 (16> 
The sought segmentation is then given by 
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$ = avg ruin ^T^ Nj • log (aj) 



(17) 



For fitting polynomials of degree two, it may be switched to a matrix 
notation. With the fitting error 



fin = s n — a/71 2 — bj7i — cj = &v. — [n 2 , y>, 1] - [a,-, bj. c.j] T 
for region Rj may be obtained 



(18) 



10 



Qj = Sj — M^-bj 
where 



(19) 



15 



g, = 



ex 
e 2 



found 







711 


1 " 














■ b i = 










1 _ 







The variance estimate is then 



(20) 



(21) 



and, minimizing &* by the estimate bj of the parameter vector, it may be 



at)., 



(22) 



20 The criterion to optimize is again given by Eq. (17), with the variance 

estimate given by Eq. (21). Extension to polynomials of higher degree is obvious. 
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For each segmentation, the optimal fit can thus be determined. The ML- 
segmentation is the one with minimum optimal fitting error. Since the number of 
possible segmentations is limited by the coherency constraint (Eq. (4)), full search is 
practically feasible. 
5 Segmentation by a Hidden-Markov Model 

The Hidden Markov model (HMM) describes a random state sequence, 
which can only be observed through another random process conditioned on it. In the 
above case, the observable random process is the 95-percentile computed from the 
vessel map histograms. The underlying, hidden process is modeled as a two-state 
10 sequence, which indicates whether or not the coronary vessel tree in an angiography 

frame is fully filled by contrast agent. It should be noted that sequences with a plurality 
of states, e.g. 3 or more, may also be modeled. The HMM is then parameterized by 

• the 2x2 state transition probability matrix A = [ay] 

• the probability density functions (pdf) pk(s), k = 1,2, of the observations 
15 conditioned on the states, and 

• the initial state probability vector II = [7C|, 712] 

Thus, the optimum state sequence Q subject to the MAP criterion is 
determined, i.e. such that P(Q|S) is maximum. To estimate the conditional pdfs pi(s) 
and P2(s), it is assumed that the histogram of the feature curve consists of two Gaussian 
20 densities. This histogram is thresholded by Otsu's algorithm such that the separability 
between the two classes is maximized. An example histogram is shown in Fig. 6. From 
this initial histogram, the threshold is found to 57. Once hard thresholding is carried out, 
mean and variance of the conditional pdfs as well as prior state probabilities 71 1 and n 2 
can be estimated. 

25 The a priori knowledge is now predominantly expressed via the state 

transition matrix A, with the transition probabilities such that staying in the current state 
is strongly encouraged. Empirically, A was determined as 



[ 0.7 0.3 1 
A I 0.3 0.7 J 



(23) 
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Two different optimization algorithms were implemented: the first 
approach uses the Viterbi algorithm to optimize over all possible state sequences, 
without considering the coherency constraint. The probability for assigning a certain 
sequence of labels Q - (q l9 q 2 ...,?*) to a certain sequence of frames l x J 2 ... 9 I N is 

with qe {/,/}, I denoting inflow and / denoting filled state. p(q ( ) is the 
probability of frame /, being of the state q t as given by the pdf that are assigning 
probabilities to percentile values (coming from Otsu's method applied to the percentile 
histograms). 
1 0 The transition probabilities 

[0.3 q t *q$-i 
are given by the transition probability matrix A . 

The computations are on the order of 4N. The coherency constraint is 
15 brought to bear afterwards by finding the state sequence consistent with the coherency 
constraint which is closest in Hamming-distance to the solution found by the Viterbi 
algorithm. The second option proceeds similar as done for the above Maximum- 
Likelihood approach by carrying out a full search over all possible sequences consistent 
with the coherency constraint. Here, the N 2 /2 different sequences need to be tested. In 
20 practice, both optimization methods gave very similar results. Note, however, that the 
Viterbi-based approach is more flexible when the coherency constraint is violated, e.g. 
by giving several bursts of contrast agent. 

Advantageously, the above described method and image processing 
device allow for an identification of complete vessel tree frames in coronary angiograms 
25 in an automatic manner. Furthermore, the above described method and operation may 
be applied to other applications, where objects, preferably moving objects of interest are 
imaged by a time series of images where the object of interest is visible due to a contrast 
medium. 



